Inter— molecular structure factors of macromolecules in solution: integral equation 

results 
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The inter-molecular structure of semidilute polymer solutions is studied theoretically. The low 
density limit of a generalized Ornstein-Zernicke integral equation approach to polymeric liquids is 
considered. Scaling laws for the dilute-to-semidilute crossover of random phase (RPA) like structure 
are derived for the inter-molecular structure factor on large distances when inter-molecular excluded 
volume is incorporated at the microscopic level. This leads to a non-linear equation for the excluded 
volume interaction parameter. For macromolecular size-mass scaling exponents, v, above a spatial- 
dimension dependent value, v c — 2/d, mean field like density scaling is recovered, but for v < v c 
the density scaling becomes non-trivial in agreement with field theoretic results and justifying 
phenomenological extensions of RPA. The structure of the polymer mesh in semidilute solutions is 
discussed in detail and comparisons with large scale Monte Carlo simulations are added. Finally a 
new possibility to determine the correction to scaling exponent W12 is suggested. 

PACS numbers: 61.25.Hq, 61. 12. Ex, 61.20.Ja 



I. INTRODUCTION 



Whereas the conformational statistics of a single flex- 
ible polymer chain in dilute and semidilute solutions are 
understood rather well, less is known about the inter- 
molecular packing. It is well understood that a semidi- 
lute polymeric solution builds up a temporary mesh with 
a mesh size, the density screening length £ p , which for 
macromolecular solutions can become large compared to 
the length scales characterizing the individual monomers 
jl], |2|. However, the inter-molecular packing inside the 
mesh but still on length scales large compared to the 
chemistry dependent local length scales is as yet unclear. 
It has been the focus of recent neutron scattering experi- 
ments |3|,H], scaling considerations and field theoretic cal- 
culations [0,§|, and of computer simulations ||. Older 
theories for the inter-molecular structure, which either 
used the random phase approximation (RPA) Q| or the 
assumption of Gaussian inter-molecular correlations Q , 
failed to incorporate the non-mean field like correlations 
on scale £ p of semidilute polymer solutions. The recent 
field theoretic results lead to contradicting results as will 
be pointed out and resolved in this contribution. 

Integral equation theories for simple liquids directly 
address the problem of inter-particle packing in dense flu- 
ids. Starting with the work of Schweizer and Curro 
this approach has successfully been extended to macro- 
molecular liquids. The polymer reference interaction site 
model (PRISM) integral equations have been fruitfully 
applied to describe inter alia the inter-molecule correla- 
tions in dense homopolymer systems, polymer blends and 
block copolymer melts |l(J . PRISM is a macromolecular 
generalization of the RISM theory of small molecules of 
Chandler and Andersen JTl , 12 . The low density limit 
of PRISM theory shall be worked out in detail in this 



contribution in order to discuss the density correlations 
on the mesh size length scale. This works either extends 



]13| [La], or complements ]16|,|17| previous studies. 

It is a priori not related to nor required for the success 
of the PRISM approach to polymer melts whether it also 
correctly captures the long ranged correlations of semidi- 
lute polymer solutions. As liquid correlations in melts 
generally are short ranged, an approach like PRISM ap- 
propriate for dense systems, need not be a useful ap- 
proach to (semi-) dilute solutions, where long ranged cor- 
relations are of interest. Nevertheless, the simplification 
of the PRISM equations to low polymer densities worked 
out here will be argued to provide a useful description 
of the inter-molecular correlations building up the poly- 
mer mesh in polymer solutions [f| |l4| . Criteria for 
the quality of the approach will be established from com- 
parisons with simulations, field theory and mean field 
results. 

The aspect of screening of the intra-molecular ex- 
cluded volume shall be neglected in this work. It would 
require the use of self-consistent PRISM theory, which is 
considerably more demanding fl(i[| . Moreover, the errors 
made, when neglecting the crossover to Gaussian intra- 
molecular correlations on length scales large compared to 
the density screening length £ p , will not affect the scalings 
of the inter-molecular correlations for distances smaller 
than £ p , which are the main focus of this contribution. 
Thus, in the following the intra-molecular correlations 
shall be characterized by a density independent polymer 
structure factor uj q , which, for macromolecules of N seg- 
ments at the positions r a , is defined as follows: 
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Its full functional form will not be required. Knowledge 
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of its variation for small, large and intermediate wave 
vectors suffices (l],||. 

For small wave vectors, the number of scattering units, 
the index of polymerization N, where N 3> f for macro- 
molecules, and the global molecular size, the radius of 
gyration R g , can be obtained from a scattering experi- 
ment measuring tv q : 



N 



N{l-- d q 2 R 



for qR g < 1 , (2) 



where d is the spatial dimension. In an intermediate wave 
vector range, the macromolecule is supposed to be self 
similar, which leads to a power law behavior in oj q deter- 
mined by the fractal dimension, dp = 1/v. 
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for l/Rg «9«1/(T 



(3) 



The fractal exponent v also determines the size-mass 
scaling, R g oc <jN v , where a smooth crossover from Eq. 
(||) to Eq. (||) is assumed around qR g « 1. The assump- 
tion of an intermediate self-similar molecular structure 
rules out the study of compact macromolecules, e. g. 
hard sphere like colloids, but is appropriate for polymer 
chains in good, v — 0.588 . . ., or 9-solvents, v = i, or 
for rods, v = 1, which share some properties with semi- 
flexible polymer molecules like actin or DNA Jl], |2], p"8| . 
The Kuhn'sche-segment size a in Eq. (|J) is of the or- 
der of local polymer-specific length scales where micro- 
scopic segmental packing effects influence the compli- 
cated structure of ui q . This chemistry dependent varia- 
tion of Lo q around qa ~ 1 can be included in PRISM stud- 
ies [[To), but shall be neglected here. Only the self scat- 
tering contribution, a = p in Eq. ([l]), which is the only 
remaining contribution for large wave vectors, qa 3> 1, is 
universal and needs to be considered. 



for qa 3> 1 



(4) 



Thus, a generic smooth crossover from the point particle 



self scattering term, Eq. (W) 
molecular correlations, Eq. (' 



, to the self similar intra- 
), will be assumed. Chem- 



istry dependent local packing will show up in all correla- 
tion functions on microscopic length scales but will not, 
except for in prefactors, affect the inter-molecular struc- 
ture on global length scales like the molecule size, R g , or 
the mesh width, as will be shown explicitly. 

In order to characterize the total, including the inter- 
molecular, density correlations of an interacting polymer 
system, further correlation functions need to be intro- 
duced. In order to compare them with results from other 
approaches it is useful to recall their definition as used 
in PRISM theory [jl0|,[l{|. To be specific, let us consider 
n polymers with TV scattering units in a d-dimensional 
volume V, where in the thermodynamic limit the num- 
ber density of segments, g — 4^, is kept fixed; i. c. 
n, V — > oo with g — const.. The local, fluctuating den- 
sity is 



e (r,t)=^^*(r-rW(i)) 



(5) 



= 1 a=l 



with equilibrium average {g{v,t)) = g. The spatial com- 
ponents of the density fluctuations shall be denoted by 
£> q , where: 



n N 



0q 



= f d d re^ g{v) = Y^Y j e 



P (i) 



(6) 



= 1 Q=l 



Their statistical average vanishes except for zero wave 
vector, (p q ) = gS^o- 

The total structure factor, S q , as measured for ex- 
ample by coherent neutron scattering [^],[l9| is given by 
the second moment of the wave vector dependent density 
fluctuations, from Eq. (0): 



{OaQq) ~ nNS^o 



nN^ q 

n N 



_L E E<e ^»-^) ) _ n ^ q0 



(7) 



ij — l a, {3=1 



The total density fluctuations are straightforwardly sep- 
arated into density fluctuations on the identical polymer, 
Lu q , the intra-molccular structure factor, and on different 
polymers, h q , the inter-molecular structure factor. 



S q — Lu q + gh q 



(8) 



where the intra-molecular part has already been defined 
in Eq. (0) above. The inter-molecular structure factor, 
h q , describes the packing of different molecules and is 
given by the restricted sum i =/= j: 



V 



N 



E E^^-^o. (9) 



The inter-molecular structure factor is the Fourier trans- 
form of the inter-molecular pair correlation function, g(r) 



(10) 



h q = / dVe iqr U(r)-l) 



The pair correlation function describes the probability 
averaged over all segments of finding at a distance r from 
site a on molecule i another segment /3 of a different 
molecule j. From Eqs. (P,|l0|), one obtains: 



lir) 



V 



N 

E E 



(*(r-(rW -■#>))). (11) 



g(r) is non-negative and approaches unity for large sep- 
arations r, because then statistical correlations between 
the sites at and have vanished |l9| . 
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Implicit in the Eqs. ([?]) to ([Tl]) is a neglect of a special 
site dependence of the density fluctuations as it might 
for example arise from chain-end effects for linear poly- 
mers (9). Star polymers, where sites in the core region 
possibly experience very different local density fluctua- 
tions than sites in the star-arms, also would require a 
more elaborate treatment jlO]]. Nevertheless, for arbi- 
trary macromolecular architectures, the above defined 
correlation functions are experimentally measurable at 
least in principle and can also be determined from com- 
puter simulations. They contain information about the 
local liquid structure and remain meaningful in the whole 
accessible density range, from dilute solutions to melts. 



II. PRISM INTEGRAL EQUATIONS 

Whereas the effects of the excluded volume on the 
intra-molecular structure, like swelling, are already taken 
into account in Eq. (Q), PRISM theory ||[h} explic- 
itly enforces inter-molecular excluded volume by requir- 
ing the pair correlation function to vanish for distances 
smaller than the segment size R. 



g(r) = 



for r < R 



(12) 



Building upon the Ornstein-Zernicke approach so suc- 
cessful for simple liquids |l9| , an averaged molecular site- 
site Ornstein-Zernicke like equation JO], [l2| is formulated 
and can be viewed as definition of an effective potential, 



the direct correlation function c q : 



h q = LO q C q 



(uj q + gh q ) 



(13) 



Equations (||) and (|l3|) can also be brought into a RPA- 
like form which supports the interpretation of the direct 
correlation function as an effective potential. 



S q =UJ q 



QCq 



(14) 



A. Thread limit for (semi—) dilute solutions 



In the simplified non-self-consistent PRISM approach, 
the intra-molecular structure ui q , from Eq. (|l|), is as- 
sumed to be given, and (mostly numerical) techniques to 
solve the integral equations, Eqs. ( ^2|) to (15) are em- 
ployed [0. Note that for a Pade approximation to to q 
for Gaussian chains with v — ^ in d — 3 an analytic so- 
lution of the PRISM equations on all length scales exists 
|]l6| . The solution technique employing the Wiener-Hopf 
factorization as pioneered by Baxter [^0| can straightfor- 
wardly be extended to Gaussian chains in (low) odd di- 
mensions, d = 5, 7, . . ., but a simpler approach can also be 
used in order to study the low density results of PRISM 
theory analytically. For the mentioned case, v — ^ in 
d = 3, this was first used in |L3],[l4|], explicitly justified in 
Jr|, and without proof extended to exponents v within 
the bounds \ — v < 2 m iF^l - Here, general arguments 
on the solutions of Eqs. (|l2[) to ( |l5| ) allow to find the low 
density limits in the more general case \ < v and d > 2. 

The special low density limit, called "thread PRISM" 
model |ll| Q , studied in the following assumes that poly- 
mer solutions can be modeled as a low density limit of 
the one-component PRISM equations for polymer melts. 
Special solvent effects, are assumed to be taken into ac- 
count via the model for the intra-molecular structure, 
Eqs. g) to (|). 

In general, the excluded volume constraint, Eq. (E3), 
and the PY closure, Eq. (|l5|), lead to a discontinuity in 
the pair correlation function g(r) and in the direct cor- 
relation function c(r) at contact: 



g d := g(r \ R) > 
c(R-) := c{r / R) ^ 



(16) 



Note, that the actual values of gd and c(R-) will depen- 
dent on details of the monomer chemistry, as the seg- 
ment size R obviously is a microscopic length. To con- 
nect both quantities it is useful to express h q and c q as 
one-dimensional Fourier transforms, 



Different from the RPA approach, c q is not considered 
to be given, but needs to be found from a solution of 
the non-linear integral equations. Besides Eqs. (|l2|) to 
(|l4|), a further equation, the "closure" approximation, 
is required to determine the solution. The most simple 
and yet appropriate closure to treat the inter-molecular 
steric repulsion is the Percus-Yevick (PY) approxima- 
tion which expresses the expectation that the effective 
potential is short ranged: 



c(r) = 



for r > R 



(15) 



Note, that the interaction described by c(r) is localized 
on microscopic length scales. Thus, the PRISM equations 
describe the interplay of local inter-molecular steric in- 
teractions and long ranged intra-molecular correlations 
due to macromolecular connectivity |p|,[Tc|]. 



h q = J dre^ ](r) , 
c q = J dr e lqr i(r) , 

with the symmetric functions, 

j(r) = n d ^i^dss(g(s)-l)(s 2 -r^ , 
i(r) - Q(R-\r\)n d ^f*dssc(s) ( s 2 -r 2 ) £ 



(17) 



(18) 



where = 



Z7T ' 



T ^ d /2) denotes the surface of the d- 
dimensional unit sphere. Because of Eqs. ( |l2]Jl5| ), i{r) 
can be non-smooth for \r\ < R only, while this can hap- 
pen for j(r) for \r\ > R. Thus, using the large wave 
vector limits, h q oc —gdCosqR/q d for q — > 00 and simi- 
larly for c q , and the large wave vector asymptote of uj q , 
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see Eq. (0), Eq. (|T^) connects the discontinuities of g{r 
and c(r) at r = R, leading to: 



9d 



c(R- 



(19) 



Moreover, one concludes that c(r) is finite. Thus, the 
scaling of the Fourier transform of the direct correlation 
function with global parameters (to be defined below) 
also can be connected to the contact value: 



c q = R d c(R~)f c (qR) =: -g d R d f c (qR) 
where the regular function, 



(20) 



f c (x) = J y<1 d d y e lxy c(Ry)/c(R), has a finite value at 
x = 0. For the analytically known results of the PRISM 
equations these property could be shown explicitly fl6| , 
and they can now be used to simplify the PRISM equa- 
tion for low densities. 

In order to extract the large distance solution of the 
PRISM equations, it proves useful to shift the micro- 
scopic length scales to zero: 



#-►0, 







with a / R = const. 



(21) 



In order to evade the trivial limit of a non-interacting 
ideal gas, the length scale of the intra-molecular correla- 
tions, to be denoted by £ c , which is proportional to the 
molecular size, £ c oc R g , is kept finite by increasing the 
index of polymerization N . 

N -> oo so that £ c cx oN v = const. . (22) 

Also, in order to keep inter-molecular excluded volume 
active, the bare segmental density is increased beyond 
bounds: 



Q 



oo so that q/ Q* = const. 



(23) 



In the thread PRISM equations, there enters a typical 
density, the dilute to semidilute crossover density 
which is familiar from scaling considerations [jjj. As will 
be shown below, g*, is defined differently for scaling ex- 
ponents v below and above a value v c , which denotes the 
crossover to mean-field like behavior. 



from Eq. p0| ) for all x. However, a solution to Eqs. 
(Ul) to ( |l5|) depending on / c (0) can be constructed on 
large distances, r R,a, or, equivalently in the thread 
limit, for finite distances, r > 0. This holds, because in 
the limit of R — > 0, the excluded volume condition Eq. 
(|l2|) affects a point (of measure zero) only, and thus h q 
can be obtained from the Fourier integral of g(r) outside 
the core, r > 0. From the inverse transformation and 
Eqs. (|l^j2^), the contact value follows in the general 
case, where R, denotes a vector of length R+: 



NA 

9d-l = 

Q 



UqR)u>, 



(2tt) 



l + Af c {qR)uj q 



(26) 



The normalized functions Co q = uj q /N and f c {x) = 
fc{x)/f c {Q) have been introduced, and the long wave- 
length interaction parameter A, which abbreviates the 
zero wave vector limit of the direct correlation function: 



.4 



-Ngc q=0 = Ngg d R d f c (0) . 



(27) 



For (semi-) dilute solutions, the limit Eq. ( pl| ) simplifies 
the PRISM integral equations, because the core condi- 
tion, Eq. (|l2]), becomes irrelevant, and Eq. ( p6] ) as- 
sures that the connection E q. ( |l9| ) is satisfied. The two 
further conditions, Eqs. (22,23), assure that non-trivial 



solutions describing an interacting polymer solution are 
obtained. Because of Eq. (^2|), only the long ranged scal- 
ing form of the intra-molecular structure factor enters: 



UJ a 



l + 0{x 2 ) x 
l/x 1 /" x 



0, 

oo 



(28) 



And Eq. ( |23[ ) enforces the molecules to interact, A ^ 0, 
so that Eq. (^) leads to a transcendental equation de- 
termining A (equivalently gj), which is the only unknown 
parameter in the thread inter-molecular structure factor: 



gh q 



-NA- 



1 + AU! 



(29) 



oc 



N 



for v <v c 



(24) 



Equations (|21j) to (|24j) specify the thread PRISM limit. 
Note, that the divergent number density actually 
corresponds to a vanishing polymer volume fraction,^), 
and thus (semi-) dilute polymer solutions are studied as 
claimed. 



g*a d oc 



]/jyOd-l) 
1/N 



for v < v c , 
for v > v c , 



(25) 



where 4> = 0(1) corresponds to polymer melts. 

The solution of the non-linear PRISM integral equa- 
tions in the general case of course demands to find f c (x) 



III. THREAD LIMIT RESULTS 

In the dilute to semidilute concentration region, the 
thread PRISM result for h q assumes the RPA like form, 
Eq. (^), where the interaction parameter A needs to be 
found from Eq. ( p6| ) . The total structure factor then also 
has simple RPA-like form: 



S a 



N- 



1 + Auj 



(30) 



where generally uj q differs from the Gaussian form and 
the integrated inter-molecular interaction strength, A, 
in general differs from the simple RPA approximation, 
A oc NgR d . Because of the large wave vector behavior 
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of Q, Eq. (28), the limit R — > affects the integral in 
Eq. (26) differently for v < v c = 1 or v > v c . In the 
first case, the integral converges uniformly, and integra- 
tion and limit can be interchanged, thus only / c (0) = 1 
enters. In the second case, the integral converges only 
because of the wave vector dependence of f c (qR), and 
leads to RPA or mean-field behavior. 



A. Below the mean field crossover 

In the thread equation for the interaction parameter 
A, Eq. d26|), the limit R — ► can be performed trivially 
for v < v c = §, and A becomes a function of g/g* only, 
where the crossover density g* — Nild/(2^ c ) d enters. 



(1 - 9d) 



g/g* 
A 



dx- 



1 + Alj(x) 



(31) 



For size-mass scaling exponents v corresponding to frac- 
tal dimensions, dp = l/v, equal to or exceeding the spa- 
tial dimension, i.e. for v < 4, the intra-molecular struc- 
ture on long length scales determines the thread equation 
@. There, screening of the intra-molecular excluded 
volume has been neglected in the present approach, and 
the polymer segregation effect predicted in the thread 
limit requires use of the self-consistent PRISM approach 
p5| . In order to avoid this complication, the exponent v 
will be restricted in the following, v > 4. From the two 



limits of the integral in Eq. (|3l|), constant for A -C 1 
and A~^ 2 ~ ud ^ for 4> 1, the scaling form for the thread 
parameter can be determined. 



Ja{x) oc 



A = (g/g*)fA(g/g*) , where 
const. + 0(x) for x — > 
for x — > oc 



2-1. 

X" d 



(32) 



In Eq. (p2]), the contact value correction was already 
neglected, as it is of higher order, as can be deduced 
from Eq. Actually, a scaling law follows from Eqs. 

(^?],[32]) for the contact value: 



(33) 



where the identical scaling function Ja from Eq. ( |32] ) 
enters. The numerical prefactor c of course depends on 
the microscopic details of the polymer model. Note that 
for v < v c the contact value vanishes like N~^~ ud > for 
N — > oo in the dilute case. The scaling function Ja also 
determines the density dependence of the mesh size or 
density screening length, as for large reduced densities, 
g 3> Q*, the width of the total structure factor can be 
estimated from S q = {N/A)/(l + (q£, P ) 1/l/ ), for q£_ c > 1, 
which leads to: 



For a given model of the intra-molecular structure 
factor, the thread equation, Eq. (|l]) with g c i = 
0, allows to determine Ja straightforwardly. Figure 
^ shows the result for the polymer model: uj q — 

J2a 0=1 ex P — 3 ~§~\ a ~ P\ 2v i with the exponent given 
by the Flory approximation v = 3/5 corresponding to 
good polymer solutions [jl|,|§]. Numerical solutions of the 
microscopic PRISM equations (|l^) to ([15]) for this model 
and for not too large degrees of polymerization, N, still 
exhibit rather large corrections to the thread asymptote. 
This can be expected to be model dependent. From Fig. 
|l| one notices that the connection of the contact value to 
the small wave vector interaction parameter, Eq. fl27|), 
already holds for values of N where the asymptotic /a, 
Eq. ([?2|) , is not yet reached. Differences appear for larger 
packing fractions and signal concentrated or melt like 
polymer packing. 

The pair correlation function, g(r), in the thread limit 
can be obtained for finite segment distances from the 
Fourier transform of Eq. (|29|). For dilute densities, 
g <C g* and thus A <C 1, it differs from the ideal 
gas limit, g(r) — 1, because of two molecule interac- 
tions. In the semidilute concentration region, g 3> g* 
and i> 1, the replacement A = (Cc/^p) 1 ^ shows that 
g(r) depends on the two length scales, £ c and £ p , indepen- 
dently. On length scales large compared to the mesh size, 
r»^, the inter-molecular structure factor exhibits the 
well known corr elation hole @,|,[l(],[l|, which asymp- 
totically for q ^> g* exactly cancels off the long ranged 
intra-molecular correlations: 



N _ 

h q ^> u(q£ c ) 

Q 



for q < l/£„ ; g > g* 



(35) 



This result is equivalent to S q <C uj q 



for q£ p <C 1 in the 

semidilute range |l]]. Note that the self-similar struc- 
ture of the molecule leads to the power law behavior 
gh q oc — (aq)~ 1 /' / for l/£ c <C q <C l/£ p , which is equiva- 
lent to a power law variation in the pair correlation func- 
tion: g(r) - 1 cx {-l/Q){a/r) d - x / v for £ p < r < £ e p). 
For larger distances, r ^> £ c , Eq. ( |35| ) describes how gjr) 
decays exponentially to its random mixing value unity. 

Within the mesh size, i. e. for distances around 
and smaller than the density screening length, the inter- 
molecular correlations do not depend on the molecular 
size and ui in Eq. (|2^) c an be replaced by its large g£ c 
asymptote from Eq. (|28|). This leads to 



-H d P 



for l/f c <C q ; g > Q* , 
(36) 



C P = c'a{ga d )-' j/{ ' yd - 1) 



(34) 



where the limiting behaviors of h q for q^ p large or small 
compared to unity can be read of immediately and 
c = N/(gA^ d ) approaches a number (c — > 6.26... for 
v = 3/5). Note, that Eq. (|34j) for the density screening 
length ensures a smooth crossover of Eq. ( |36"| ) to Eq. ( |35| ) 
in the correlation hole region. The variation of the pair 
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correlation function, which describes the mesh structure 
for r <C £c: can thus be obtained in closed form, if the 
neglect of the cutoff of the correlation hole at r 15 £ c in- 
cluded in Eq. (j35[), is kept in mind. For r <C £ C) g( r ) 
depends on r/£ p only, with: 



g(r) = 1 - c 



d d y 



-iyr/(,p 



4. v 2/v 



(2ir) d 

1 - {i P /r) d - 1/v r » £ P ; r « & , 



(37) 



numerical constant c, which however may depend on the 
ratio a/R of the microscopic length scales. 



A = cg/g* 



(39) 



In the semidilute density regime, the width of the total 
structure factor again determines the density screening 
length, S q = (N/A)/(l + (q£ P ) 1/u ) for q£ c > 1, with the 
result: 



H P = c' a(ga d ) 



(40) 



where c ensures 17(0) = in agreement with Eq. ( |3l| ) 
and constant prefactors of order unity have been sup- 
pressed in the final two lines. For polymer chains in good 
solvents, the smooth increase, g(r <C £ P ) oc (r/^p) 1 ^ 3 , 
by accident agrees with the estimate from Ref. [ MJ , 
g{r -C £ p ) oc (r/Cp)^~ 1 ^ 1 ' ~ (f/Cp) 1 ^ 3 where 7 is associ- 
ated with the entropy of a single polymer chain Q . The 
depth of the correlation hole displays an intriguing de- 
pendence on the fractal and spatial dimensionalities. The 
probability to find a segment of another polymer within 
the considered molecule decreases strongly if 1/v — > d. 
From Eq. @ one estimates g(r w f c ) - 1 oc -1/N ud ~ x , 
which becomes a number of order unity in the case 
v = 1/d. The smooth variation of g(r) at short distances 
explains, why the scaling of the correct contact value gd, 
Eq. ([33j), with macroscopic variables can be estimated 
from the thread solution, Eq. (37), by gd oc <?(<x); its 
dependence on the ratio of the microscopic length scales, 
c r/ ' R, however cannot generally be recovered in this way 
jig . Note that Eqs. ( p4| ) to (|3^) asymptotically apply for 
semidilute solutions, g 3> £>*, whereas ( p9| , |30| , |32"| ) describe 
the full dilute-to-semidilute crossover region. 



B. The mean field cases 

The condition Eq. ( p6| ) for the contact value gd, or 
equivalently, for the thread parameter A becomes inde- 
pendent of the microscopic interaction details only for 
v < u c . Above the crossover exponent, v > v, 



= 7v the 



integral over the effective potential as it enters h q is de- 
termined by the local structure in f c (qR). In the thread 
limit Eq. (|2^) becomes a linear, density-independent 
equation for gd with solution: 



9d = l/[l + 



d d q 



{2TT/R) d 



e i ^{qa)- 2 l»f c {qR)\ 



(38) 



Thus, a finite density independent contact value follows 
in the mean field like cases v > v c . Obviously, its ex- 
act value depends on the solution of the PRISM equa- 
tions considering all microscopic details and is beyond 
the reach of the thread PRISM approach. The interac- 
tion parameter A thus shows the density scaling as ex- 
pected within RPA. The reduced density g/g* appears, 
with g* = l/(Na d ), and A becomes — with a unknown 



Again, the numerical prefactor, c', depends on the poly- 
mer model. The result, Eq. (|35|), discussed for the inter- 
molecular structure in the correlation hole region holds. 
The mesh structure factor, h(q ^> l/£ c ), however, shows 
a different density scaling: 



-co d {tplofl" 
{qip) 1/v + {q£,p) 2/l - 



for 1/& « q , (41) 



which again, with Eq. ([io|), leads to a smooth crossover 
for intermediate distances, l/£ c -C q <C 1/£ P , but to a 
small distance divergence of the thread pair correlation 
function, g(r) oc —{<j/r)( d ~ 2 /^ for r -C £ p . This results 
from the neglect of the wave vector variation of the direct 
correlation function, Eq. (^) , and reinforces that the va- 
lidity of the thread g (r) is restricted to r 3> cr for v > v Cl 
where g(r) is still positive, as it must be by definition. 



IV. DISCUSSION AND COMPARISON WITH 
OTHER APPROACHES 

In this work, scaling limits appropriate for the dilute 
to semidilute concentration regime of macromolecular so- 
lutions have been derived starting from the microscopic 
PRISM integral equations. RPA like expressions for the 
total density fluctuations, the structure factor S q , Eq. 
d30| ) , were given, where the density scaling of the interac- 
tion parameter, A, was deduced from the local excluded 
volume constraint, Eq. (12). The thread interaction pa- 
rameter is connected to the more familiar excluded vol- 



ume parameter, v, via A 



_ N \M g 

- -ur 



v, where Na is Avo- 



gadro's number, M the molecular and Mo the monomer 
weight. Effective density dependent excluded volume pa- 
rameters v(g) have often been used in connection with 
RPA expressions ||pi]], and Eq. (32) justifies this. 

The crossover of the single chain correlations to Gaus- 
sian large distance behavior for r 3> £ p due to intra- 
molecular excluded volume has been neglected and would 
affect the model for uj q , Eq. (||), and consequently the 
thread results for large distances. Use of self-consistent 
PRISM Q to incorporate this would be required, but 
Eq. (32) for the thread parameter A indicates that no 



change of its density scaling can be expected. 

The crossover density g*, Eq. (|24|), arises from the 
full microscopic PRISM equations as the relevant low 



G 



density scale, and importantly, the qualitatively differ- 
ent definitions in the mean field, v > v c = §, and in 
the non-trivial cases, v < i/ c , are recovered. Whereas 
for v < v c the molecular crossover density, c* = g*/N, 
is defined in terms of the molecular size only, c* oc 1/ R d 
for v > v Cl in the mean field cases also a microscopic 
length, the segmental hard core diameter R, enters, 
c* oc l/(R 2 g /u R d ~ 2 ^) for v > v c . This indicates that 
inter-molecular steric interactions become important as 
soon as the macromolecules fill space for v < v c , whereas 
for the more open molecules, v > v Cl much higher densi- 
ties are required. 

For chain polymers, the upper critical dimension, 
which separates mean field and fluctuation dominated 
structures, agrees with the renormalization group results, 
d c = 2/v c = 4 jjjgj. For rod polymers, v = 1 > v c , the 
mean field like behavior underlies the successful Onsager 
theory of the nematic transition j2^] and is generally ar- 
gued to be true f||]. Note that in the studied PRISM the- 
ory orientational, "nematic" , interactions are not treated 
correctly pof , and thus a nematic transition for rods is 
missed. Very recently, PRISM has been generalized to 
treat oriented polymer fluids and the isotropic-nematic 
liquid crystal transition |24j . As the (isotropic) crossover 
density c* for rods is of the order of the nematic transition 
density J22|, a suppression of nematic order is required 
to study the described isotropic semidilute rod solutions 
experimentally. Networks of stiff semiflexible molecules 
like actin may provide good systems |lS| ] . 

Of course, the full PRISM integral equations, which 
have been introduced to study dense polymer systems 
with short ranged melt-like correlations |9], incorporate 
wave vector dependent corrections in e. g. the effective 
interaction c q , see Eq. flTi|), when l/q approaches local 
length scales. 

From the compressibility, which is connected to the 
zero wave vector limit of the total structure factor, the 
equation of state can be obtained, where II denotes the 
osmotic pressure |lQ,[L9|: 



n 



i 



gksT 



Q Jo N 
±(l + 0(g/g*)) g< e . , 

d 



(42) 



QCT 



when A is given by Eq. ( |32] ) . The non-mean field behav- 
ior for v < v c for semidilute concentrations |l5t agrees 
with the exact Des Cloizeaux result and the second 
virial coefficient, 1/{Nq*) oc R d g /N 2 , recovers the picture 
of dilute polymer coils interacting like hard spheres of ra- 
dius R g but it does not vanish for dilute O-solvents, 
i. e. for v = ^ in the present approach. PRISM the- 
ory apparently correctly captures the leading asymptotic 
behaviors, the free molecule limit II = (g/N)kBT for 
g <C g* and the TV-independent power law for g ^> 
but the next to leading terms are not described correctly 
in general. 



The structure of the polymer mesh in semidilute so- 
lutions, i. e. the inter-molecular structure factor, h q 
, on length scales of the order of the density screening 
length, £ p , has not been conclusively discussed from first 
principles calculations. The thread PRISM results for 
the non-mean field like case of polymer chains in good 
solvents give explicit results, Eqs. (g^,^), which can be 
compared to results from other approaches. 



A. Comparison with scaling considerations 



Detailed scaling law considerations of h q in the limit 
qR g 3> 1 have been presented in Q and can be di- 
rectly compared with Eq. (|36|). The limit, h q oc for 
1/Rg -C q <C l/£ p strongly differs from the correlation 
hole behavior, gh q = —(qcr) 1 ^ with a the Kuhn'sche seg- 
ment size, predicted by the thread PRISM theory for this 
wave vector window. Physically, the long ranged varia- 
tion of h q arises from the rearrangement of the polymer 
mesh around a molecule on distances up to the molecule's 
size. This adjustment compensates for the excess den- 
sity due to the considered molecule, leading to the small 
total density fluctuations expected for concentrated sys- 
tems. The correlation hole has first been predicted and 
discussed for polymer melts but PRISM theory also 
predicts it for semidilute solutions || , in agreement with 
scaling considerations in j2l]] but in disagreement with 
the mentioned scaling picture presented in p|. Intrigu- 
ingly, PRISM theory, Eq. (|37j), recovers the tendency of 
macromolecules to segregate for v = 1/d as discussed for 
ideal chains in two dimensions |5J|. 

In agreement with the thread PRISM result, a scaling 
law is postulated in Ref. pi] for the intermolecular struc- 
ture factor inside the coil radius, h(q ^ 1/R g ) = h(q^„ 
which leads to the prediction h q oc q~ d for q — » oo 



As a scaling law can only hold for distances large com- 
pared to the microscopic length scales, « 1, this 
result can be compared with the thread scaling power 



law, h n 



oc q 



in Eq. (BO), and again differs. 



Within thread PRISM the behavior of h q arises natu- 
rally as it matches smoothly to the microscopic limit, 
h q ~ H cos qR for qR ^> 1, because the contact value 

vanishes asymptotically, oc (a 7£ P ) 2 / Iy ~ d for £ p 3> a. 
This supports the expectation in EjJ . Computer simula- 
tions could address this question for polymer chain solu- 
tions as shown in section IV. C, where corrections to the 
low density scaling law, Eq. (p3|), need to be considered 
which will arise due to finite packing fractions. 

Accepting the existence of a scaling law for the contact 
value of macromolecules in solutions, then Eq. ( ^3|) can 
be used to connect the PRISM results to field theoretic 
calculations for two-polymer systems. 
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B. Comparison with field theoretic calculations 

Field theoretic calculations which employ the mapping 
of the polymer problem onto the 0(n — > 0) magnetic 
model lead to numerous single chain results and have re- 
cently been extended to provide information about the 
inter-molecular structure factor on short length scales 
In Q, the mentioned behavior h q oc q~ d for q — > oo 
is recovered from the field theoretic calculation and used 
to support the scaling picture discussed in the previous 
section. The implications for g(r) can be compared with 
another field theoretic calculation which studies the num- 
ber of intersections of two polymer chains Let £ 2 (i? e ) 
be the number of intersections of two random (or self- 
avoiding) walks whose end-to-end distance is R e : 



^Mm= E<*(tf>-tf>-B.)*(r« 



V(4ira 2 ) d / 



))• 



a, 0=1 



(43) 



For intermediate distances R ei a -C R e <C R g , where the 
two polymers overlap but local effects do not dominate 
£2, the scaling £2 oc {a/R g ) UJl ^ p > is predicted, where 
the two molecule correction to scaling exponent Wi 2 ap- 
pears ||. The contact value can now be obtained from 
£2 by integrating over all possible distances and (triv- 
ial) factors of normalization, as can be seen from Eqs. 

©nil). 



9d 



d d R e 



£ 2 (-R e 



(Aira 2 ) d / 2 N 2 



(44) 



Using the results for two polymers from M to obtain the 
scaling of the contact value in the dilute case, one finds: 



T?d-U} 12 (P) 

„RG fll 

9d oc N2 



for g — > 



(45) 



The thread PRISM result, Eq. (||), differs from this in 
general, because in PRISM theory the correction to scal- 
ing exponent is approximated to w*2 rcad = 0. Its value in 
quadratic order in e — 4 — d is known, and the value ap- 
propriate for polymer chains in good solvents turns out 
to ui 12 (G) = \e- ±fe 2 + . . . ~ 0.40 @, which can be 
compared to the PRISM and to the mean field approxi- 
mation, ix>^> PA = d — 2/v, which qualitatively differs be- 
cause it is negative. The thread approximation, u>\2 = 0, 
is correct at and above the upper critical dimension d. 

It appears difficult to envisage simple forms of g(r 
which reconcile the prediction h q oc q for q — > 00 
with the results for £2 Q , and the dilute limit of thread 
PRISM theory qualitatively agrees with the later. 



i 



C. Comparison with Monte Carlo simulations 

Monte Carlo simulations are well suited to study the 
inter-molecular structure of polymer solutions but face 



the difficult challenge to achieve a clear separation of the 
three length scales, segment size a (or excluded volume 
size R), density screening length £ p , and molecular size 
R g (or molecular correlation length £ c ); see the discus- 
sion in H . Whereas in fi in the range 1 / R g <C q <C l/£ p 
a discrimination of the two predictions, h q oc const, from 
scaling considerations [Q and h q oc q~ x with x w 1/v 
(as follows from the PRISM treatment of the correla- 
tion hole) appears possible and appears to support the 
latter, no clear conclusions about the exponent in the 
asymptotic behavior, h q — » q~ x for l/£ p <C q -C l/R, 
with x = d = 3 (scaling picture), x = 4 (RPA), or 
x = 2/u « 3.34 (thread PRISM) were possible. Even 
recent large scale simulations of the bond fluctuation 
model (BFM) [^5],|2(| do not provide a conclusive test 
of the large q dependence if h q is considered |^7j] . Figure 
U shows data from Ref. J27|] for semidilute solutions and 
rather large chain lengths, N = 2048, where £ c = 94, 
£p = 14 for c/c* = 97.7 and £ p = 31.1 for c/c* = 25.6, 
and the steric segment size is R = 2 in units of the lat- 
tice constant of the BFM. For a fit, the asymptotic thread 
PRISM prediction, Eq. (|3"6|), is shifted by a factor indi- 
cating that non-asymptotic corrections to c cannot be 
neglected. 



A clearer picture of the polymer mesh structure is 
provided by the pair correlation function g(r), which is 
the Fourier transform of h q and asymptotically should 
follow Eq. (|37| ) in the thread limit, a,R <C £ p and 
i P *C ic oc R g . Figure || shows the two g(r) for the above 
parameters, where £ p is defined by collapsing the simula- 
tion data onto the master curve at g{r = £ p ) = 0.747. 
Note that this is an unfamiliar definition of £ p which 
gives values (theoretically) proportional to the standard 
ones. These values of £ p also produce the collapse of the 
h q onto a common curve shown in the inset of Fig. ^ 
and lead to a reasonable collapse of the pair correlation 
functions onto a common master curve. Finite size cor- 
rections enter from short distances because of the finite 
excluded volume segment sizes, R/£ p . These corrections 
can also be understood as finite packing fraction correc- 
tions. Large distance deviations from a common curve 
appear because of the finite chain sizes, £ c /£p- Never- 
theless, the short distance behavior of g(r) provides a 
sensitive test of the various predictions. The prediction 
h q —> l/q d would correspond to a logarithmic variation 
of g(r — > 0), which appears to be ruled out by the data. 
Also the thread PRISM prediction, g(r -> 0) oc r 1 / 3 , 
appears incompatible with the data, even if finite seg- 
ment size corrections are approximated by a shift of the 
r-origin. The RPA prediction for Gaussian polymers, 
g(r — > 0) — gd oc r, of a linear increase in r, can describe 
the data over small intervals (like 0.1 < g < 0.25) but 
fails to account for the slight curvature of especially the 
lower density curve. Moreover, the contact value of the 
RPA cannot be expected to vanish asymptotically if pa- 
rameters appropriate for a fit to h q are used. An increase 
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in the range of a power law fit to g(r) at the lower den- 
sity up to an interval 0.1 < g < 0.42 is possible if the 
following assumption about the pair correlation function 
for semidilute solutions is made: 

9(r) - g(r/£ p ) for £ p -> oo ; £ p /£ c - , (46) 
g(x < 1) oc x 2/»-d+"i2(P) . (47) 

where for dilute cases the same power law with the re- 
placement £ p — > £ c can be expected from scaling consid- 
erations. This power law would match the scaling law 
for g(r) for r — > smoothly to the calculated vanish- 
ing contact value gd from Eq. (fl5|). Note that such a 
matching is predicted by PRISM. In Eq. ( p7|) however, 
the exponent is corrected because the correction to scal- 
ing exponent w\2(P) is taken into account. According to 
scaling arguments ||[28|, there is a term of the form of 
Eq. ( |47| ) present in the intra-molecular correlations also, 
although it is masked by chain-end effects there. 

The expected power law for good solutions, g cx r om 
for r < £p with v = 0.588 and w = 0.40 is com- 

patible with the simulation data, if a finite shift owing 
to a finite segment size is anticipated. The power law 
h q — » l/g 2 / i/ + Wl2 ( p ) also is compatible with the data as 
shown in Fig. ^, but could less be argued on data for h q 
only. 



V. CONCLUSIONS 



The thread PRISM results derived and discussed here 
justify earlier phenomenological extensions of RPA like 
expressions. The density dependence of the excluded 
volume parameter is derived from a microscopic incor- 
poration of inter-molecular excluded volume and intra- 
molecular connectivity. Various comparisons with rigor- 
ous field theoretic calculations show that leading asymp- 
totic predictions, even for non-mean field like situations, 
are captured correctly in the PRISM integral approach. 
The correction to scaling exponent, which appears in 
the molecular mass dependence of the contact value of 
two polymers, provides a typical example where thread 
PRISM provides a much better description than mean 
field theory but fails to describe all non-trivial correla- 
tions. PRISM theory suggests useful concepts like the 
pair correlation function g(r) and predicts scalings laws 
which provide a framework for the interpretation of data 
if the exponents are corrected. Thread PRISM thus turns 
out rather useful for semidilute solutions, where it ex- 
plicitly describes the inter-molecular correlations of the 
polymer mesh and results from more rigorous approaches 
are scarce. Moreover, as PRISM theory is successful for 
polymer melts, it provides the unique possibility to ap- 
proach polymer systems at all densities. 
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FIG. 1. 

Scaling function /^(c/c*) (bold solid curve) of the thread interaction parameter A versus rescaled (molecular) 
polymer concentration in double logarithmic presentation for the Flory exponent, 1/v = 1.67; c* = l/(2ir 2 Q) is the 
molecular overlap concentration. The curves labeled with degree of polymerization, N, are results from microscopic 
PRISM calculations of the model described in the text, where £ c = 0.28N u a is found with (psj). Full symbols give 
/a determined from S q —>o, and open symbols from the contact value gd shifted by model dependent factors (0.088, 
0.086, 0.085, 0.085 with increasing TV). 



FIG. 2. 

Inter-molecular structure factors h q from a Monte Carlo simulation of the BFM for polymers of length N = 2048 
(£ c = 94 for v = 0.588; all lengths given in units of the lattice constant of the BFM). The data are taken from j27j and 
are shown scaled with £ p determined from the corresponding g(r)'s of Fig. ||. The +'s belong to the rescaled density 
c/c* = 25.6 and the x's to c/c, = 97.7. The asymptotic thread PRISM result, Eq. ([56]), shifted by a correction factor 
0.25 is shown as solid line, whereas the dashed line indicates a power law, q -3 - 80 following from Eq. (jrif). Bold solid 
lines mark where q = 1. 



FIG. 3. 

Pair correlation functions g(r) versus rescaled distance r/£ p of the BFM for the two densities, c/c» = 25.6 (thin 
solid line, £ p = 31.1) and c/c» = 97.7 (thin dashed line, £ p = 14.0) from |27|. The choice of £ p collapses the curves 
at g = 0.747. Circles mark the contact values gd- The asymptotic thread PRISM prediction, Eqn. ( |3^ ) (thick 
solid line), and shifted according to a finite segment size (thin dot-dashed curve) are shown. A small r asymptote, 
g = 0.99((r — 0.855i?)/£ p ) ' 80 according to Eq. ( ]47j ) is indicated by a long dashed curve, where R — 2 in units of the 
lattice constant is the excluded volume segment diameter of the BFM, and a linear asymptote is given by a dotted 
line. The bold solid line denotes £c/£p for the lower density. The inset enlarges the small-r region showing the shifted 
thread PRISM, the power law and the linear curves with the same line types as in the main part. 
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Fuchs and Miiller, Inter-molecular ... Fig. 1 




Fuchs and Miiller, Inter-molecular ... Fig. 2 




Fuchs and Miiller, Inter-molecular ... Fig. 3 



